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ABSTRACT 


The East Greenland Current and its associated Marginal 
Ice Zone is a region of intense dynamical activity. A two 
layer, primitive equation, numerical model is used to 
simulate an eddy-jet interaction in the East Greenland 
marginal ice zone region. The effects of wind direction, 
topography, and sense of eddy rotation on the eddy-jet 
interaction are examined to determine the seaward ice 
transport, icebanding, and dipole formation. It. Xs 
determined that an anticyclone (15 cm/s) interacting with 
a jet (30 cm/s) will develop a dipole that advects ice away 
from the ice edge. The dipole formation and ice advection 
away from the ice edge is not seen for a cyclone-jet 
interaction. It is also seen that a jet with no winds 
flowing parallel to the ice edge will create an iceband due 
to the cross ice edge Ekman transport. The interaction of 
both the cyclone and anticyclone with the jet creates 
downstream perturbations in the jet leading to a sinuous 
ice edge. Winds greater than 10 m/s dominate the ice 


dynamics over that induced by the ocean flow fields. 
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I. INTRODUCTION 


A. THE EAST GREENLAND CURRENT REGION 

The East Greenland Current and its associated marginal 
ice zone (MIZ) is a region of intense dynamical activity 
with complex mesoscale ocean processes influencing the ice 
edge. Figure 1.1, from Paquette et al. (1985) shows the 
general circulation of the region. Since the Fram Strait 
region is the primary source for Atlantic Water flowing 
into the Arctic basin and Polar water flowing out of the 
Arctic Basin, the region's importance both tactically and 
Sementirtically is striking. Experiments such as the 
Marginal Ice Zone Experiment (MIZEX), have been conducted 
this decade using shipboard sensors, Lagrangian drifters, 
and satellite imagery. From the results of these studies, 
an understanding of the mesoscale dynamics of the 
oceanographic phenomena is beginning to evolve. Figure 
1.2 is a remarkable photograph of a cyclonic eddy/ice edge 
interaction in the marginal ice zone. The marginal ice 
zone is divided into two different regions of ice 


concentration, Gascard et al. (1988). The outermost region, 


along the ice ocean interface, is delineated by a sinuous 
edge with an ice concentration of 3/10 or less. The inner 
region is defined by an ice concentration of 34/10 or 
greater.  Gascard et al. (1988) summarizes the MIZ as an 
area of low ice concentration, typical floe sizes ranging 
from one to 25 km, and variable eddy features located 


between these two limits. 


B. EDDY CHARACTERISTICS 

A number of mechanisms, five of which were proposed by 
Johannessen et al. (1987), may be responsible for the 
generation of the observed mesoscale eddies. Among the 
generation mechanisms are barotropic and baroclinic 
instability of the East Greenland Current or marginal ice 
edge front, differential Ekman pumping along the ice edge, 
topographic generation, and ice edge instability driven by 
internal ice dynamics. Johannessen et al. (1983) 
documented a 10 km wide ice edge mesoscale ocean front in 
the East Greenland Current. The jet front was located in 
4000 meters of water with velocities of 10 cm/s. Several 
mesoscale eddies with length scales ranging from 5 to 15 
km were also documented. Both the jet and the eddies 


exhibited little vertical shear. Based on observations 


and theory, Johannessen et al. (1983) concluded that the 


eddies resulted from barotropic instability of the oceanic 


meont. The energy derived by the eddies came from the 
kinetic energy of the mean flow. Pagure 1.39 from 
Johannessen et al. (1987) depicts schematically the 


interaction of an open ocean mesoscale eddy with the ice 
edge. 

According to Gascard et al. (1988), from observations 
during the 1983 and 1984 marginal ice zone experiment, the 
East Greenland Current (EGC) is a narrow, stable current. 
Gascard et al. (1988) suggested that the West Spitzbergen 
Current (WSC) and the Norwegian current are the main 
generators of mesoscale eddies in the Fram Strait. The 
eddies are then advected by the return Atlantic flow 
following  bathymetric fracture zones into the East 
Greenland Current. OC RS SE 2 (1938) Support this 
argument by concluding the lateral turbulent heat flux 
associated with the East Greenland Current is insufficient 
to be consistent with local generation of eddies through 
@roelinie instability. 

Using Lagrangian drifters, Gascard et al. (1988) traced 


this eddy path and documented the interaction of open ocean 


eddies with the ice edge and the East Greenland ocean front 
(Figure 1.4). While some of the eddies were anticyclonic, 
the majority were cyclonic, consistent with the earlier 
findings of Johannessen (1983). Because the East Greenland 
Current follows closely the continental slope of East 
Greenland, topography may also be important in the 
interaction of the eddies with the East Greenland Current 
and the MIZ. 

Wind forcing has been shown to cause upwelling and 
downwelling at the ice edge. Hakkinen (1986b) showed how 
winds parallel to the ice edge with the ice edge on the 
right would produce upwelling. Conversely, a wind with 


the ice edge to the left would produce downwelling in the 


Pyucnoelimes wind at the ice edge can also result in 
vorticity input into the ocean. This is demonstrate Sn 
Smith et al. (1988) where upwelling favorable winds were 


shown to destroy ocean anticyclones at the ice edge. Along 
ice edge winds may thus also be important in the 
interaction of eddies with the East Greenland Current 
marginal ice zone. These processes will be discussed in 


further detail in the ne t sece len: 


C. PURPOSE 

"uuEthrs study, the interaction of an open ocean eddy 
with a marginal ice zone and associated along ice edge 
ocean front is considered. Results are obtained uSing a 
coupled air-ice-ocean numerical model. The effects of 
winds and topography on the eddy interaction with the MIZ 


and the East Greenland ocean front are included. 
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Figure 1.2 MIZ Ice Edge Cyclonic Eddy 


Source: Johannessen et al. (1987) 
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Figure 1.3 Schematic Ice Edge- Eddy Interaction 
Observed During MIZEX 84 


Source:  Johannessen et al. (1987) 











Les 


* 


Surface Drifter Trajector 
Eddy Trajectories across 
(Dots represent 0000 UT on each 


Lagrangian 
Indicating 


trait 
day) 


1.4 


Figure 


Fram 


Gascard et al. (1988) 


* 
* 


Source 


2 


II. PREVIOUS WORK 


Numerous mesoscale physical processes have been 
identified as important in the MIZ region. A brief 
description Of previous findings will facilitate 


explanation of results in this study. 


A. UPWELLING/DOWNWELLING 

Winds varying with time have distinct dynamical effects 
on the ice edge. Hakkinen (1986b) showed how winds 
parallel to the ice edge, with the ice edge on the right, 
would produce upwelling. The momentum flux of the air to 
the ocean is smaller than the momentum flux from the air 
to ice to ocean so continuity dictates upwelling in the 
pycnocline MEE gore 2.15 The opposite holds true for a 
wind with the ice edge to the left and there would be a 
downwelling band. With the assumption of a sinuous ice 
edge, Hakkinen (1986b) showed how Ekman pumping drives 
both upwelling and downwellina along the ice edge, Fiaure 
2.2. It is also demonstrated how eddies are formed at the 
ice edge through differential Ekman pumping. This study 


will show how a sinuous ice edge can occur from the 


10 


weer action Of a cyclonic or anticyclonic eddy and a 


Stable barotropic jet. 


B. OPEN OCEAN EDDY INTERACTION WITH THE MIZ 

Smith et al. (1988) showed that at high latitudes an 
eddy decays much more slowly than at lower latitudes. 
This is because the westward propagation speed of the eddy 
is limited to the Rossby long wave speed, c = - RB R,?, 
where c is the wave speed, f is the planetary vorticity 
gradient, and R, is the Rossby radius of deformation. At 
80°N, the Rossby long wave speed (c) is on the order of 
.01 km/day, so the decay time associated with the Rossby 
wave radiation takes much longer. Since the simulations 
in this study are limited to periods of less than six 
days, the observed eddy motion and decay are not due to 
Rossby wave dynamics. 

The interaction of a mesoscale open ocean eddy with a 
marginal ice zone in the absence of an along ice edge 
ocean front is considered in Smith et al. (1988). They 
determined that the Coriolis term in the ice momentum 
equations, without ocean surface slope terms, caused 
divergence of the ice over a cyclone and convergence of 


the ice over an anticyclone. When the pressure force in 


vi 


the ice caused by the ocean slope is included in the ice 
momentum equations, a geostrophic balance occurs in the 
ice. This means that the pressure gradient force is as 
important as the Coriolis force and should be included in 
the equations of motion. The geostrophic balance holds 
true for a light wind situation only. When winds are 
greater than 10 m/s, the geostrophic balance becomes 
secondary and the wind forcing dominates. Smith et al. 
(1988) also demonstrated that upwelling winds provide 
cyclonic vorticity to the ocean destroying anticyclopme 
eddies (see Figure 2.3). SAME EN yy downwelling winds 
provide anticyclonic vorticity which destroys cyclos 
eddies. They also suggested that since 12 of the 14 
eddies observed by Johannessen et al. (1987) were cyclonic 
eddies, the upwelling winds during MIZEX-84, destroyed the 
anticyclonic eddies. In this study, both upwelling ama 


downwelling winds will be examined. 
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Figure 2.1 Ekman Induced Upwelling (W,) by Along Ice 
Edge Winds (symbols defined in Appendix) 
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Figure 2.2 Wind Forced Upwelling/Downwelling in a 
Sinuous Ice Edge 


Source: Hakkinen (1986b) 
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mure 2.3. Cyclonic Vorticity Input to Upper Ocean by 
Along Ice Edge Winds (AC denotes existing 
Ocean anticyclone) 
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C. EDDY JET INTERACTION 

The interaction of an eddy with an ocean jet in the 
absence of ice was first examined by Stern and Flierl 
(1986). In their study, the interaction of a point vortex 
with a jump discontinuity jet was considered on an f- 
plane. Their experiments illustrated a new eddy 
propagation tendency not associated with planetary 
rotation or background advection. If an anticyclonic eddy 
approaches normal to a jet, the negative vorticity of the 
eddy interacts with the cyclonic, positive vorticity side 
of the jet, and a weaker cyclonic eddy is formed. The 
propagation tendency of the eddy pair is determined by the 
addition of the two vorticities and the dipole is advected 
in the common direction. 

This interaction of an eddy with a Jet resulting in 
eddy motion was also demonstrated in Smith and Davis 
(1989), Figure 2.4, where eddies and jets with finite 
vorticity distributions were also considered. The extent 
to which these processes apply in a marginal ice zone 
Situation is examined in this study. 

In the previous eddy-jet interaction studies of Stern 


and Flierl (1986) and Smith and Davis (1989), the motion 


16 


Of a vortex away from a Jet was found for eddies initially 


within a nondimensional length scale, R(o), from the edge 
of the jet. The nondimensional length scale depends on 
the vorticity of the eddy and its size. Smith and Davis 


(1989) found eddy motions for values of R(o) in the range 
BENE tO 2.5, with insignificant eddy motion for greater 
R(o). values. For comparison with Smith and Davis (1989), 
eo), in this study, is approximately 1.1 indicating that 


vortex motion induced by the jet is likely to occur. 
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Figure 2.4 Eddy-Jet Interaction 


in the Absence of Ice (a) 
Upper Layer Relative Vorticity (b) Lower Layer 


Relative Vorticity (C.I = .05x10"s” for Day 
0-6 and C.I.: = »02x210. = Zforpasal 2 BS 


Source: Smith and Davis (1989) 
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D. ICEBANDING 

Icebanding is a frequently observed dynamic phenomenon 
that occurs in the marginal ice zone. Icebanding is 
important because it causes a more rapid ablation of the 
ice margin according to Hakkinen (1986a). Hakkinen 
(1986a) examined several possible mechanisms for iceband 
nation and their typical size. First indications from 
satellite images are that the icebands are 10 km wide 
with approximately 10 km of ocean separating them. Later 
measurements taken by surface platforms indicate narrower 
icebands between .5 to 1 km. One of the theories cited by 
Hakkinen (1986a) is that internal waves cause icebanding. 
The wavelengths and phase speeds of internal waves in a 
two layer system are similar to the motion of the 
icebands. It is also suggested by Hakkinen (1986a) that 
previous studies demonstrate the formation of icebands by 
lee waves. 

Hakkinen (1986a) also demonstrated the formation of 
icebands by modeling the ice edge with ice-ocean being 
coupled through interfacial stresses. The momentum flux 
used by Hakkinen (1986a) has an air-ice coupling three 


time greater than the air-ocean coupling. Along ice 
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(upwelling) winds were introduced that varied sinusoidally 
with a four day period. The along ice edge velocity would 
drive the cross ice edge Ekman flow and some icebanmdia 
would develop in 12 days. To have ice banding occur at 
least one wind reversal was necessary. This study 
illustrates similar results in a much shorter time by the 


introduction of an ocean jet along the ice edge. 


20 


III. NUMERICAL TECHNIQUE AND MODEL PARAMETERS 


The numerical model was initialized with parameters 
chosen to resemble the East Greenland Current with a 
velocity of 30 cm/sec as observed by Gascard et al. (1988). 
A schematic of the model domain is presented in Figure 3.1. 
An ice edge in which concentration increased to the west 
was included and both cyclonic and anticyclonic eddies of 
10 km diameter were examined in a time sequence.  Foldvik 
et al. (1988) determined a typical tangential current speed 
in the observed eddies to be 15 cm/s . The initial eddy- 
jet interaction experiments are without winds. The effects 
of along ice upwelling and downwelling winds on the eddy- 
jet interaction are also considered. Topography simulating 
the continental slope off Greenland was also included in 
Some runs to determine the effect on the eddy jet 
interaction. The determination of the various parameters 


will be discussed in chapter four. 


zu 


A. THE OCEAN MODEL 

The model used in this study is the same as that used 
by Smith et al. (1988) with the addition of a barotropic 
jet. Experiments are performed using a two-layer primitive 
equation semi-implicit numerical scheme. In “mer 
experiments the simulated ocean topography has a flat 
bottom. In several experiments the topography of the East 
Greenland Slope is represented by a linear bottom slope, 
sloping upward toward the west. The depth ranges from 4050 
meters to 1050 meters at the boundary. The effects ot 
topographic dispersion of the eddy are examined. 

The model constants C, C „r Cw, etc. are chosen from 
Hakkinen (1986b) and are defined in the Appendix. The 


motion of each layer is governed by the momentum equation: 


aV, z^ e =, n. ô aw w 
a ya A OS + 417”) (eqn. 3.1) 
and the continuity equation: 
ol (eqn. 3.2) 


"an uM PV, = 0 
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for layer (i-1 upper and i-2 lower) thickness h,, and 


velocities v,. The fluid is hydrostatic and Boussinesq. 
Fluid density (Pi) in each layer is fixed and no mixing 
across the fluid interface is allowed. The ocean-ice 


coupling is through the ice-water interfacial stress, 


ri" 


= Baer tg — thy) | 1 — thy | 


(eqn. 3.3) 


and concurrently, the ocean-air coupling through the air- 


water stress, 1°": 


a == MEN, E Uy) | Ug 77 My (eqn us. ) 


for ice (uj), air (u,) and upper layer ocean (u,) velocity 
vectors. The small-scale turbulent eddy dissipational 
processes are represented by a horizontal  Laplacian 
operator on transport. The coefficient À, has been chosen 
small (10m?s") making the fluid relatively inviscid, since 
the Laplacian operator only roughly represents a complex 
process. All notation is defined in appendix. 

The semi-implicit scheme has been used in a number of 


mid-latitude ocean mesoscale circulation studies. A 
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thorough explanation by Hurlburt and Thompson (1980) and 
Smith and O'Brien (1983) provide more details of the 
scheme. They show how the model conserves mass and total 
energy in the absence of dissipation. The model 
verification of the Rossby dispersion characteristics was 
done by comparing linear analytic solutions to the linear 
test cases run on the model as discussed in Smith and Reid 
(1982). The characteristics of the upwelling and 
downwelling response at the ice edge are consistent with 
that seen in the models of Hakkinen (1986b), and Smedstad 


and Roed (1985). 


B. BOUNDARY AND INITIAL CONDITIONS 

A rectangular (111 x 79 km; i-111,  J=/79) en 
difference gridded domain is used for the model. In order 
to align the domain with the East Greenland Current, it is 
rotated approximately 90 degrees relative to the 
referenced latitude of 80? N. This simulates a southward 
jet running approximately parallel to the east coast of 
Greenland. The initial state consists of a 
geostrophically balanced jet extending uniformly across 
the basin from north to south and a Gaussian eddy in 


gradient balance east of the jet. Experiments are 
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Pitetalized with Dbareeropie.  velocmey dissaibutlons. 
While the East Greenland Current velocities at depth may 
be weaker than near the surface, this choice eliminates 
unstable meandering of the jet, simplifying the eddy-jet 
anmseraction. This also allows a more direct comparison 
with the previous eddy-jet study of Smith and Davis 

LOREN 
The barotropic eddy and jet are defined in terms of 
upper (h,) and lower (h.) layer thickness. 
= (ean 3.5) 


and in the lower layer: 





en uM" 
hy) = Thy — Ate P=), ze ) 
2m au 


(eqn. 3.6) 


where a negative value is assigned for cyclones and a 
positive value for anticyclones. H, and H, are mean layer 
thickness values. 


The eddy radius is defined as follows: 


reS E a) ES a (eqns: 3T) 
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The eddy center is located at the x,, y, coordinates. The 
eastern edge of the jet is at y, (y grid point j = 38), 38 
km from the eastern boundary. In this study the initial 
eddy location relative to the edge (y=y.) of the jet, is 
fixed at ten km to the east. XE. and L, are the c- foi mmi 
width scales (= 5 km) for the jet and eddy. The amplitude 
of the Gaussian jet was chosen to give a maximum velocity 


(v of approximately 30 cm/s in each layer in the base 


aoa) 
case experiments. This velocity is similar to observed 
velocities of the East Greenland Current (Gascard et al. 
1988). A radiation condition was used on the downstream 
(left) boundary in all model simulations. The radiation 
condition, (Camerlengo and O'Brien 1980), advects flow out 
of the basin at speed 3 when flow is outward adjacent 
to the boundary. On the upstream boundary (right), 
constant jet inflow of the form of the initial condition 
is continually maintained. The north and south boundaries 
are no-slip walls where both tangential and normal flow 
are set equal to zero. The initial condition parameters 
are provided in the appendix. The thermocline depth in 


the East Greenland Current ranges from 25-100 meters. The 


upper layer mean thickness was chosen to be 50 meters in 
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all cases. The lower layer mean thickness is 4000 meters. 
This gives a first internal Rossby radius of deformation (R,) 
equal to 5 to 7 km, consistent with the observations of 
Johannessen et al. (1983). Nondimensional eddy size y= 


L methus order 1. The Rossby number for this flow ( vw... ) 


IA X 


is .21 so the nonlinear terms can not be neglected in the 


momentum equations. 

In the wind driven simulations, spatially uniform winds 
are specified (10m/s) along the ice edge. As in previous 
studies, Hakkinen (1986b), Smith et al. (1988), the wind 
vector has been rotated to offset cross ice edge drift 
associated with the ocean Ekman drift to the right. The 
angle of rotation which gave the minimal meridional ice edge 
drift was found to be 25° to the left of ice edge 
orientation. The simulations are integrated for duration of 
6.5 days, within which significant eddy jet interactions 
occur. The eddies and jet are initialized with barotropic 
structure but baroclinic structure in the interface is free 


to develop during the model run. 
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ee THE ICE MODEL 

The ice model is initialized so the ice edge runs 
parallel to the jet. The motion of the ice is governed by 
the momentum equations (eqn. 3.8 and 3.9) and the 
continuity equations (eqn.3.10 and 3.11). Although the 
ice is initially at rest, adjustments to an equilibrium 


state occurs rapidly, normally one-half day. 


Og On, Ou iiti O(/t, + hy) 

NOE Tue DR Al ud dum ER = fv, + — ONE Bp vec (eqn. 3g) 
Ov, Ovi CUN ai iw O(h + I) 
Ay) eq" 

A O(mu)  O(mv,) ) 
Op ucc E (eqn. 3.10) 
Qu. ULM al en We 

Ot Ox Oy (eqn. 3.11) 


The ice concentration is represented by (A) and the ice 
mass per unit area by (m). The air-ice coupling for the x 
component is represented using equation 3.12 and for the y 


component using equation 3.13. 


ai 
mm psu. = uj) | Ua — uj 


(eqn. 3.12 
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AUTRES 5 
Cy ENE Vy) | Va — v] (eqn E 15) 


Ice thickness (D) is computed using equation 3.14, with an 
initial thickness of 2 meters. The concentration is 
initialized at 05$ at the open ocean boundary and 75$ at 


the western boundary, Figure 3.2. 





(0141) (eqn. 3.14) 


Most mesoscale studies neglect the pressure gradient 
force in ice associated with the sea surface slope, but 
Smith et al. (1988) showed that it should be included in 
the ice equations in no wind or light wind studies. The 
Laplacian damping term has also been included in the 
continuity equation (eqn. 3.10), as in previous studies 
such as Hakkinen (1986b) and Hibler (1979). Since dynamic 
effects are the focus of this study, no internal stress 
terms are considered, an assumption which is appropriate 


mamice concentrations less than 85% (Hibler 1979, 1984). 
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Figure 3.1 Coupled Ice Model and Initial Conditions 
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Figure 3.2 Ekman Divergence in Ice Concentration Over 
a Laterally Sheared Ocean Jet 
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D. WIND FORCING 

In some experiments wind forcing is included to 
determine the effect on the eddy-jet interaction. The 
winds are either upwelling favorable (ice edge to the 
right) or downwelling favorable (ice edge to the left). 
The predominant wind direction is seasonally dependent but 
during MIZEX-84, Johannessen (1987) observed predominantly 


northerly and southerly flow, Figure 4.1. 
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IV. EXPERIMENTS 


A. PURPOSE OF EXPERIMENTS 

The purpose of the following experiments is to 
understand the effect of eddy-jet interactions along an 
ice edge. Observed dynamic features such as 
upwelling/downwelling, dipole formation, icebanding, and 
ice transport (ice tongues) will be examined using the 
numerical model. The output of the numerical model runs 
produced a time sequence, .5 day interval, of the ice 
concentration, upper level relative vorticity, lower level 
relative  vorticity, surface height anomaly and the 


interface height anomaly. 


B. PRELIMINARY EXPERIMENTS 

A list of primary parameters is provided in the 
appendix. Table 1 gives a summary of the parameters varied 
during each experiment. 

Before examining the eddy-jet interaction in the MIZ a 
base case experiment is conducted to examine the effects 


of a jet running parallel to and under the ice edge. The 
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eddy was not included in the model initial conditions for 
this run. 

1. Simulations with No Wind 

a. Experiment 1A  (Jet/Ice edge simulation) 

In this simulation a 30 cm/s jet was 
initialized under and parallel to an ice edge without 
winds. The experiment is limited to a duration of 6.5 days 
by computing constraints, however this is sufficient to 
examine the eddy-jet interaction processes of this study. 
It is important to note that icebanding has only been 
Simulated in numerical modeling by a time reversing wind, 
Hakkinen (1986a) or upwelling favorable winds , (30° off 
ice edge with ice to the right) Smedstad and Roed (1985). 
This simulation illustrates that iceband formation can 
occur without winds. 

A time sequence of the ice concentration field for the 
base case, day 1, 2, and 4 is shown in Figure 4.2. By day 
four a decrease in ice concentration above the jet can be 
observed. In order to quantify these data a cross section 
normal to the ice at midbasin of ocean u velocity and ice 
concentration was taken for each day. On day 1, the ice 


concentration over the jet shows a decrease from .75 to 
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.60 over a 5 km band in the domain, (Figure 4.3). The 
minimum ice concentration is located over the maximum u 
merocity of 30 cm/s. On day 2 (Figure 4.4), the ice 
concentration over the jet shows a decrease in ice 
concentration from greater than .80 to .40 and then back 
to .80 over the 10 km maximum u velocity area. This 50% 
decrease shows substantial banding beginning to occur in 
just two days. By day 4, (Figure 4.5), the ice 
concentration shows a decrease from greater than .80 to 
less than .40 over the u maximum velocity 10 km area. It 
is also evident that the ice band is beginning to move 
slightly westward of the jet. 

Since there is no wind in this simulation, the cause of 
the ice banding is the ocean induced Ekman transport in 
the ice. Since the u maximum of a laterally sheared jet 
has a maximum Ekman transport (to the right in the 
northern hemisphere), a divergence occurs over this 
maximum, Figure 3.2. This effect tends to transport ice 
normal to the u velocity maximum creating an iceband to 
the left. The formation of an iceband is a function of 
the ocean u velocity and the larger the ocean u velocity 


maximum, the greater the ice v component or Ekman 
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transport. It appears that by day 4 the entire iceband is 
beginning to drift with the Ekman transport across the 
domain. 

Although Smedstad and Roed (1985) were able to simulate 
an iceband by upwelling favorable winds in 12 days, this 
study is different in that the numerical model includes a 
surface tilt term in the momentum equations to balance the 
Coriolis term. The inclusion of the surface tilt term 
results in a smaller iceband than that seen in Smedstad and 
Roed (1985). In Smedstad and Roed (1985), the along ice 
edge upwelling favorable winds create the jet which then 
causes divergence in the ice. In this study, the jet 
simulating the East Greenland Current, in the absence of 
winds, causes the divergence. 

Another interesting result is that for a jet of 15 cm/s 
or less the icebanding did not develop in 6.5 days. The 
jet has to be stronger than 15 cm/s to create sufficient 
transport to cause a divergence in the ice during this 
period. 

2. Simulation with Wind 

Smith et al. (1988) conducted a study of the effect 


of upwelling winds on ocean eddies along an ice edge but 
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did not consider the effect of a Jet running parallel to 
the ice edge. 


a. Experiment 1B (Jet ice edge simulation with 
upwelling favorable winds) 


The initial conditions in this simulation are 
the same as in the previous run; however upwelling 
favorable winds are added. Based on the findings of 
Smedstad and Roed (1985), it was expected that the 
upwelling winds would supplement the  icebanding by 
augmenting the upwelling and thus the iceband. A cross 
section of ice concentration and u velocity, Figure 4.6 
indicates that this is not the case. The concentration 
decreases less than 10% on day 4. The upwelling band is 
present at the interface by day 6, Figure 4.7, but the 
iceband is not as pronounced in the jet only case. This 
is consistent with the findings of Smith et al. (1988). 
Along ice edge 10 m/s winds contribute more to the ice 
momentum balance than does the ocean. Thus while the 
ocean jet causes some divergence in the ice, the dominant 
velocity is wind driven to the south. Computational 
considerations stopped the model at day 6 so a longer run 
time may have shown the iceband occurring as in Smedstad 


Enc Roed (1985). The important inference is that a jet, 


S 


such as the East Greenland Current may create an iceband 
on time scales as short as one day, thus contributing to a 


more rapid ablation of the ice edge. 
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TABLE 1. SUMMARY OF EXPERIMENTS 


ee nm nn en M m — ee 


EXP # JET EDDY WIND WIND TOPOGRAPHY 
(15cm/s) DIRECTION SPEED 


-—— — Ze aaa 


1A 30 cm/s none none none flat 

1B 30 cm/s none upwelling 10 m/s flat 
favorable 

2A 20 cm/s AA none none flat 

2B Ed em s -CYC. none none flat 

3A 30 cm/s A/C upwelling 10 m/s fiat 
favorable 

3B 30 cm/s A/C downwelling 10 m/s flat 
favorable 

4A MA C: upwelling Ems flat 
favorable 

4B SOR eh seme vc. downwelling 10 m/s flat 
favorable 


upwelling favorable is 205° relative to north 
downwelling favorable is 025° relative to north 
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Figure 4.1 Fram Strait, Mean Wind Direction and Speed 
(30 minute averages) 


Source: Johannessen et al. (1987) 
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Figure 4.2 Ice Concentration Days 1, 2, and 4 


(Jet-Ice Edge with No Wind, Experiment 1A) 
Consour Intervalgs05 


(Dimensions in this and subsequent domain 
contour plots are 79 x 111 km) 
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Figure 4.3 u Velocity and Ice Concentration AcrosSQ MES 
Day 1 (Solid line - ice concentration ans 
dashed line - u velocity) 

(Jet-Ice Edge with No Wind, Experiment 1A) 
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Figure 4.4 u Velocity and Ice Concentration Across MIZ, 


Day 2 Solid line Ce concentration and 
dashed line - u velocity) 


(Jet-Ice Edge with no wind Experiment 1A) 
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Figure 4.5 u Velocity and Ice Concentration Across MIZ, 
Day 4 (Solid line - ice concentration and 
dashed line - u velocity) 

(Jet-Ice Edge with no wind Experiment 14A) 
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meaqure 4.6 wu Velocity and Ice Concentration Across MIZ, 
Day 4 (Solid line 


ice concentration and 
dashed line - u velocity) 


(Jet-Ice Edge with upwelling winds, 
Experiment 1B) 
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Figure 4.7  Interfacial Height Anomaly (meters) Day 6 
(Jet-Ice Edge with Upwelling Winds, 
Experiment 1B) 
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C. PRIMARY EXPERIMENTS 
1. Experiment 2A (Base Case: Anticyclone-Jet no wind) 
The anticyclone/jet interaction in the marginal ice 
zone was chosen as the base case to demonstrate the 
production of a dipole from a monopole vortex. The 
experiment begins with the anticyclone adjacent to the 
cyclonic (and seaward) side of the laterally sheared Jet 
(see Figure 4.8, day 0). It has been shown by Smith and 
Davis (1989) that when an eddy interacts with a stable 
barotropic zonal jet, vortex pairing occurs and the eddy 
propagates away from the jet. Figure 4.8 demonstrates 
that this process can occur for East Greenland Current 
parameters as well. AS is seen in the upper layer 
pouCJEy field, the anticyclonic vorticity from the eddy 
momemeeracting with the cyclonic vorticity side of the 
jee. This interaction causes a small cyclonic vortex to 
separate from the jet. The dipole persists through day 6 
propagating away from the ice edge. This propagation is 
consistent with the interaction of opposite vorticities 
which contribute to a common propagation direction. 
The movement of the ice edge in response to the 


interaction is shown in Figure 4.9. The first indication 
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of an interaction is the perturbation created on day 1. 
The newly formed cyclonic vortex is advecting ice out away 
from the ice edge creating an ice tongue. Observations of 
the ice tongues are not uncommon (Gascard 1988). As the 
experiment continues to day 6, the ice tongue extends 
seaward approximately 10 km and is centered primarily over 
the. cyclonic vortex. This could be one of the reasons 
mostly cyclonic ice swirl eddies are observed in imagery. 
The Synthetic Aperture Radar Imagery from MIZEX 1987, 
Figure 4.10, shows the advection of ice off the ice edge 
is a manner consistent with eddy-jet dipole formation. 
Another interesting phenomena from this base case is 
the sinuous ice edge that occurs from the initial eddy-jet 
‘interaction. Hakkinen (1986b) showed that if the model 
was initialized with a sinuous ice edge, differential 
Ekman pumping would spin up anticyclonic and cyelonse 
eddies in the upper ocean. This study shows that another 
mechanism can cause a sinuous ice edge is an eddy-jet 
interaction. The initial perturbation in the Je 
advected downstream with the mean flow and another dipole 
occurs from this interaction. Figure 4.8 demonstrates how 


the perturbation in the cyclonic side of the jet causes a 
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cyclonic vortex to spin up downstream which interacts with 
the anticyclonic side of the jet causing an anticyclonic 
vortex to spin up adjacent to it. The ice concentration 
as seen in Figure 4.9 is highly variable downstream 
creating a complicated ice edge and open areas in the 
interior ice pack within six days. Ice perturbations 
downstream of eddy-jet interactions are also documented 
in the observations of Gascard et al. (1988). 

Gascard et al. (1988) observed that the initial 
perturbation caused by a cyclonic eddy interaction with 
the East Greenland Current produces two intermediate 
cyclones downstream with a separation distance of 90 km. 
Since the period of an eddy is two to three days and the 
East Greenland current has a velocity of 30 cm/s, 90 km 
is a the approximate distance a parcel will travel in that 
period of time. Gascard et al. (1988) also discussed 
eddies with a wavelength of 50 km, one-half of the 
previous wavelength. This is attributed to eddies 
developing a two ice tongue system with one ice tongue to 
the north and one to the south. The southern ice tongue 
is related to the off ice entrainment due to the cyclonic 


rotation and the northern ice tongue from a divergence 
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area upstream. The wavelength of the perturbations 
observed in this study is approximately 90 km and the 
downstream perturbations are observed for both cyclone and 
anticyclone-jet interactions. 

2. Experiment 2B (Cyclone-Jet without winds) 

Twelve of the fourteen eddies observed by 
Johannessen et al. (1987) in the MIZEX 84 data were 
cyclonic.: Thus the only parameter in this experiment 
differing from experiment 2A is the sign of the eddy 
FOLIE ome The results, however are markedly different. 
Figure 4.11 shows the vorticity distributions in the upper 
layer. On day 1 the cyclone has just started to perturb 
the jet. By day 4 the jet is slightly sinuous and by day 
6, the initial perturbation is advected to the downstream 
boundary. 

Figure 4.12 shows the effect of the cyclone-jet 
interaction on the ice concentration: On days 1 and 2 
there is a weak ice tongue being advected cyclonically 
away from the ice edge. By day 3 the Jet appears to 
advect the ice tongue downstream and smooths it out. On 
days 4, 5, and 6 the ice tongue is gone and only the 


perturbation from the cyclone can be seen. The reason ice 
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is not advected away from the ice edge as in experiment 2A 
is because the cyclonic eddy is interacting with the 
cyclonic side of the jet. Since both vorticities are of 
the same sign, no dipole is formed and no propagation 
mechanism exists to advect ice away from the ice edge. 
The eddies in experiment 2A or 2B do not have any effect 
on the icebanding; the iceband forms as it did in the 
preliminary experiments but tends to conform to the 
sinuous pattern caused by the eddy interaction. 


J. Experiment 3A (Anticyclone/Jet with Upwelling 
Winds) 


This experiment will examine the effect of 
upwelling winds on the  anticyclone-jet interaction. 
Figure 4.13 shows the dramatic effect of the wind on the 
upper level vorticity. In the anticyclone-jet with no 
winds, experiment 2A, the dipole developed rapidly, but in 
this experiment with upwelling winds the process seems to 
be damped by the winds. It is observed that a weak dipole 
develops but by day 4 the anticyclone has become very 
weak. Figure 4.14 shows the lower level vorticity 
pattern and the characteristic dipole observed in 
experiment 2A is readily evident. Figures 4.15A and 4.15B 


contrast the velocity vectors in the no wind experiment 2A 
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with this experiment which has upwelling winds. In the 
upwelling wind case, an upwelling band develops on the 
eastern boundary due to boundary conditions. This creates 
a 10-15 cm/s jet opposing the anticyclone, thus damping the 
anticyclone in comparison to the no wind case. Figure 4.16 
shows the effect of the interaction on the ice edge. The 
initial perturbation is advected downstream and by day six 
the ice edge is relatively smooth. An iceband comparable 
to the jet-only upwelling wind case, Experiment 1B is seen 
in this experiment as well. 

Figure 4.17 characterizes the interface height 
anomaly between the two layers and by day 6, the interface 
has upwelled 10 meters height in the middle of the domain. 
AS in previous studies, along ice edge winds with ice to 
the right cause this upwelling. It is also apparent in the 
figure that the upwelled interface is steeper under the 
higher concentration ice. Hakkinen (1986b) showed that 
when pycnocline changes are similar to upper layer 
thickness, non-linear effects can be important. In the 
case of upwelling winds, the Ekman driven flow advects the 
upwelling band iceward providing asymmetry. This Ti 


illustrated in Figure 4.18 by examining distance vs. time 


Sy 


it). At day 5, there is a gradual slope up to the 14 
meter pycnocline anomaly then a sharp drop as the ice edge 
is neared at x = 100. It is readily apparent that the 
slope is gradual up to the 10 meter interface height 
anomaly and steep on the ice edge side of the upwelling 
band. This phenomenon has also been observed in the East 
Greenland Current marginal ice zone by Johannessen et al. 
(1983). Hakkinen (1986a) also stated that the wind forced 
jet caused the ice to move faster at the ice edge than it 
moved inside the ice pack. 

In this experiment ice is not advected away from 
the ice edge in the presence of upwelling winds because a 
wind, greater than or equal to 10 m/s will dominate the 
ice dynamics. The cyclonic vorticity at the ice edge, 
induced by upwelling favorable winds, adds to the cyclonic 
wemereity of the jet and subtracts from the anticyclonic 
wemercity of the jet, Figure 4.19. The anticyclone is 
initially at j=30 to the south of this wind induced 
e loonie vorticity. Its decay thus appears unrelated to 
the wind. The eastern boundary develops an upwelling band 
with a jet opposing the anticyclone. The anticyclone is 


damped by the artificial jet before a strong dipole can 
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form. The initial ice perturbation is advected downstream 
with the mean flow. The upwelling wind makes the ice edge 
less sinuous. 


4. Experiment 3B (Anticyclone/Jet with Downwelling 
winds) 


An anticyclone-jet experiment is conducted to see 
the effects of a downwelling favorable winds, (ice edge to 
the left), on the eddy-jet interaction. Since the winds 
are opposing the jet, the jet is substantially damped in 
the upper layer as evidenced by the decreasing gradient in 
the surface height anomaly, Figure 4.20. The winds 
provide anticyclonic vorticity to the jet, damping gene 
cyclonic shear needed for dipole formation. There is no 
dipole formation and the characteristic ice edge of 
experiment 2A on day 2 contrasts vastly with the ice edge 
on day 2 of this experiment. Once again the wind forcing 
dominates the ice dynamics of the experiment. 


5. Experiment 4A (Cyclone-Jet Interaction 
with Upwelling favorable winds) 


This experiment demonstrates chat upwelling 
favorable winds damp the cyclone-jet interaction in the 
upper level but generate a downstream cyclonic eddy due to 


differential Ekman pumping. This experiment is first 
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mermaucted using a velocity of 30 cm/s for both the eddy 
and the jet. The stronger eddy causes a stronger 
perturbation in the ice edge hence making the ice edge 
more  sinuous, Figure 4.21. The upwelling band is 
interrupted by a downwelling band in the sinuous ice edge, 
Figure 2.2, appearing to spin up a downstream cyclone with 
a 15 meter amplitude seen in the interface, (Figure 4.22). 
Whereas experiment 2B, cyclone-jet interaction with no 
wind, shows a sinuous ice edge and upper level cyclonic 
vorticity out to day six, this experiment shows the upper 
level cyclonic vorticity weakening by day five, Figure 
4.23. The cyclonic vorticity is however stronger than the 
corresponding upwelling wind anticyclone case. The ice 
edge is less sinuous than the no wind case and the initial 
perturbation is advected downstream with the mean flow of 
the jet. This case was run again with an eddy with a 
=en onal velocity of 15 cm/s and a jet of 30 cm/s. 
Similar dynamics are observed; however, the stronger eddy 
creates the stronger perturbation which enhances the 
downstream eddy formation. 

An important aspect of this simulation is that the 


winds blowing along the sinuous ice edge create a 


TO 


downstream cyclonic eddy as evidenced by the 15 meter 
interface height anomaly. The generation of this 
downstream cyclonic eddy is consistent with the generation 
theory of Hakkinen (1986b). The eddy-jet interaction 
provides the mechanism for the sinuous ice edge assumed in 
the Hakkinen (1986a) study. Another interesting aspect of 
this mechanism for ice distortion is that unstable 
meandering of the jet is not required. 

6. Experiment 4B (Cyclone-Jet with Downwelling winds) 

This experiment is conducted to contrast with 
experiment 3B and experiment 4A. The cyclone-jet with 
downwelling winds did cause the formation of an small ice 
tongue during the initial interaction in contrast to the 
anticyclone-jet with downwelling winds (Figure 4.24). 

The primary difference between this and experiment 4A 
is that the downwelling winds oppose the jet and inhibit 
the advection of the ice tongue downstream with the mean 
flow. The ice tongue remains stationary through day two. 
The downwelling band of five meters in the interface 
height anomaly supports the expected results from the 


downwelling favorable winds discussed by Hakkinen (1986b). 


56 





DAY 3 


Upper Level Relative Vorticity (* 10°) Day 
0-6 (Anticyclone-jet with no wind, 
Experiment 2A) 


Figure 4.8 


5 





DAY 1 










Å met 
Fe une 
x > . ie 






mn == 
===" — 


= å 
== n 
i “ie? NN () . 70 f 


1. 


me — me 
+ 


== 
= = E 
x 





DAY 4 DAY 6 


Figure 4.9 Ice Concentration for Day 1, 2, 4, and 6 
(Anticyclone-jet with no wind, Experiment 
2A (Contour interva ROS) 
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Figure 4.10 Synthetic Aperture Radar Imagery MIZEX 1987 
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Figure 4.11 Upper Level Relative Vorticity (*10°) Day 
1, 2, 4 and 6 (Cyclone-Jet no winds, 
Experiment 2B) (Contour interval 5.0 * 
TO) 
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EMgure 4.12 Ice Concentration for Day 1-6 
(Cyclone-Jet no winds, Experiment 2B) 
(Concour interval- 05) 
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Figure 4.13 Upper Level Relative Vorticity, Day 1-6 
(Anticyclone-jet with upwelling winds, 
Experiment 3A) 
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DAY 6 


Eure 4.14 Relative Vorticity in the Lower Level, 
Day 6 (Anticyclone-jet with upwelling 
winds, Experiment 3A) 
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Figure 4.15A Velocity Vectors and Speed Contours (cm/s) 
at the Surface Day 6  (Anticyclone-jet no 
winds, Experiment 1A) (Contour interval 
= 2.5 cm/s) 
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Figure 4.15B Velocity Vectors and Speed Contours (cm/s) 
at the Surface Day 6  (Anticyclone-Jet 
with upwelling winds, Experiment 3A) 
Contour interval = 5 cm/s) 
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Figure 4.16 Ice Concentration for Day 1-6 
(Anticyclone-Jet with upwelling winds, 
Experiment 3A) (Contour interval .05) 
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Figure 4.17 
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Figure 4.18 
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Figure 4.19 A Relative Vorticity in the Upper Layer 
Across the MIZ, Experiment 1A 
Jet/No Eddy/No Wind Day 4 Along I-20 
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Figure 4.19B Relative Vorticity in the Upper Layer 
Across the MIZ, Experiment 2A 
Jet/No Eddy/No Wind Day 4 Along I=20 
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DAY 0.5 





DAY 2.5 


Figure 4.20 Height Anomaly in the Upper Level, Day 0.5 - 
Das (Anticyclone-jet with downwelling 
winds, Experiment 3B) (Contour inrerval ts 


meters) 
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DAY 5 DAY 6 


Figure 4.21 Ice Concentration Day 1-6 
(Cyclone-Jet with upwelling winds, 
Experiment 4A) 
(Contour Interval .05) 
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Figure 4.22 





DAY 6 


Interfacial Height Anomaly, Day 2, 
(Cyclone-Jet with upwelling winds, 
Experiment 4A) 


71 





DAY 1 DAY 2 





DAY 3 DAY 4 





DAY 5 DAY 6 


Figure 4.23 Relative Vorticity in the Upper Level, 
Day 1-6 (Cyclone-Jet with upwelling winds, 
Experiment 4A) 
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Figure 4.24 Ice Concentration Day 1-2 
(Cyclone-Jet with downwelling winds, 
Experiment 4B) (Contour interval .05) 
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7.  Topography 

Many experiments were run using the topographic 
slope of the East Greenland region to examine the effect 
of the slope on the eddy-jet interaction. Since the runs 
were all barotropic, the jet and the eddy interact with 
the bottom. Most of the effects discussed above are 
largely unmodified by the inclusion of topography. Dipole 
formation for an anticyclone-jet interactions can still 
occur, however, however, there is some evidence of 
topographic dispersion contributes to the dipole decay. 
The downstream growth of perturbations in the jet induced 
by the eddy are less over a sloping bottom. This is 
consistent with vorticity Conservati onk The topogramPr 
thus stabilizes the jet. For the purpose of this study, 
topography is seen to have no significant effect on the 


eddy-jet interaction in the upper level. 
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V. CONCLUSIONS 


The East Greenland Current and its associated marginal 
ice zone is a region of intense circulation. The purpose 
of this study was to examine the eddy-jet interaction in 
this region. The understanding of the dynamics of this 
area and their effects on the ice edge is important to 
naval operations as well as commercial shipping. 

A two layer, primitive equation, numerical model was 
used to simulate the eddy-jet interaction while varying 
selected parameters. The effects of wind direction, eddy 
rotation, and topography on the interaction were examined 
to determine the dipole formation, seaward ice transport, 
and icebanding. The model output consists of ice 
concentration, upper and lower layer relative vorticity, 
surface height anomaly, and interfacial height anomaly. 

Several icebanding mechanisms have been modeled in the 
past (Hakkinen 1986a) using wind forcing. This study 
demonstrates how an ocean jet generates an iceband by Ekman 
transport within one day, while wind forcing generally 


takes up to twelve days. The eddy interaction with the jet 


T9 


has no effect on the banding other than to cause MUS 
conform to the sinuous pattern created during the initial 
interaction. 

The development of an ocean dipole is observed in the 
anticyclone-jet interaction due to the interaction of the 
anticyclonic vorticity of the eddy interacting with the 
cwelonzeverticrity ¡OCN re A second smaller eddy 
develops from the jet and this dipole advects ice away from 
the ice edge due to the propagation tendency away from the 
ice edge. The dipole development is not observed in the 
cyclone-jet interaction case. 

The observation of ice tongue formation by Gascard et 
al. (1988) is also observed in the experiment involving 
the interaction of an anticyclone-jet with no winds. On 
day six there is a substantial ice tongue developing from 
an originally straight ice edge due to the perturbation of 
the jet by the anticyclone. The ice tongue has 
anticyclonic curvature, curving toward the downstream 
boundary. Another experiment in which there was ice tongue 
development was the cyclone-jet interaction with 
downwelling winds. Although the ice tongue is not as 


pronounced, it is interesting to note that the ice tongue 
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has cyclonic curvature toward the upstream boundary. TE 
was seen that the interaction of a eddy and a jet perturbs 
the ice edge and the perturbation was advected downstream. 
The distance between ice edge perturbations is 
approximately 80-90 km, similar to that observed by Gascard 
et al. (1988). It should be noted that the model has a 
constant jet velocity of 30 cm/s while the East Greenland 
current is a variable 30-40 cm/s, Gascard et al. (1988). 

The interaction of an eddy with a jet in the marginal ice 
zone is also important in that it is a mechanism for the 
formation of a sinuous ice edge. The sinuous ice edge is 
shown to favor wind generation of ice edge eddies through 
differential Ekman pumping (Hakkinen, 1986b). Tuo S EIS 


observed in the experiments with upwelling favorable winds. 


Winds perform an important role in the dynamics of the 
eddy-jet interaction. The inportant conclusion is that 
the winds 10 m/s or greater dominate ice dynamics in the 
upper layer. Inzbetch “the Zantrevelones and) eyelöne-jet 
interaction with winds of 10 cm/s the ice edge becomes less 
sinuous. while dipole formation can still occur in the 


case of anticyclone-jet interaction with upwelling winds, 


my 


there is no ice export seaward suggesting that dipole 
events seen in SAR imagery are most likely from light wind 
periods. The upwelling and downwelling band development 
from along ice winds is observed in the wind experiments. 
This supports the upwelling and downwelling along ice 
theory of Hakkinen (1986b). 

The same experiments conducted with a flat bottom were 
repeated using a linearly sloping bottom to a maximum depth 
4010 meters, simulating the East Greenland ocean 
topography. Topographic dipole formation and dispersion in 
the interfacial height anomaly are observed but the effect 


of topography on the eddy-jet interaction was small. 


D 


VI. RECOMMENDATIONS 


Observations have shown that the East Greenland Current 
has substantial vertical shear. Since this study evaluates 
the effects of a barotropic eddy-jet interaction in the 
marginal ice zone, it would provides a good basis for a 
follow on study for a barotropic eddy interaction with a 
baroclinic jet. 

Additional experiments with a more realistic bottom 
topography of the region should be conducted to examine 
the effects of the bottom topography in a baroclinic jet 
barotropic eddy interaction. Time dependent, off ice edge 
winds should also be included in the study to simulate 
those winds observed by Johannessen et al. (1987) (Figure 
4.1) to examine their effects on icebanding. Similar 
experiments using a larger domain to eliminate the effect 
of the artificial jet on the dipole formation should be 


examined. 
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APPENDIX 


SYMBOLS AND NOTATION 


Ice Concentration 

Eddy maximum amplitude 
Laplacian lateral friction 
coefficient for ice 


concentration 


Laplacian lateral 
friction coefficient 


Jet maximum amplitude 


Laplacian lateral friction 
coefficient for ice mass 


Variation of Coriolis parameter 
with latitude 


Air-ice interfacial stress 
coefficient 


Air-water interfacial stress 
Coefficient 


Ice-water interfacial stress 
coefficient 


Phase speed of perturbation wave 
Ice thickness distribution 
Grid spatial resolution 


Time increment 
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30 


10m?/s 


30 


3.8 x 


DD T 


m/s 
m/ p, A 
IO 


600 s 


j 0 Tm S -1 


ig 


Por 


TO 


Kronecker delta function = 


Coriolis parameter for = 
mean latitude 


Q when i-2 


id cuc Suc 


Gravitational acceleration = 9,8 m/s° 
Nondimensional eddy size IER 
Upper layer mean thickness = 50m 
Lower layer mean thickness = 4000 m 
Instantaneous layer thickness = 50 m 
e-folding scale of the ocean eddy = 5 km 
e-folding scale of the ocean jet = 5 km 

Ice mass per unit area = kg/m’ 
Pressure in the upper layer =g (Nith. td) 
Pressure in the lower layer = JE SG s 
Nondimensional eddy strength Ev Re 
upper, lower layer potential = a. 

g r ticity h; 

First internal Rossby radius 

of deformation =[g"H,H,/f," (H,+H,) 1" 
Rossby number Sv FE 
Density of air = 1 kg m" 
Density of ice EOS Eco ms 


Air-ice interfacial stress vector 


Air- water interfacial stress vector 


Ice-water interfacial stress vector 
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V 


max 


X, Y 


V 


i 


Eddy maximum tangential vector 


Cartesian coordinates directed 
N and W respectively 


Upper, lower layer relative 
vorticity 


Gradient operator 


Laplacian operator 
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Il 


Il 


15 cm/s” 
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